Quantum Simulation of Cooperative Jahn-Teller Systems with Linear Ion Crystals 
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The Jahn-Teller effect explains distortions and non-degenerate energy levels in molecular and solid-state 
physics via a coupling of effective spins to collective bosons. Here we propose and theoretically analyze the 
quantum simulation of a many-body Jahn-Teller model with linear ion crystals subjected to magnetic field 
gradients. We show that the system undergoes a quantum magnetic structural phase transition which leads to a 
reordering of particle positions and the formation of a spin-phonon quasi-condensate in mesoscopic ion chains. 
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Introduction - Quantum many-body physics is motivated 
by the description of solid-state systems, where a variety of 
complex intriguing phenomena emerge as a result of strong 
correlation effects. Our understanding of the latter is hindered 
by the intrinsic computational complexity of many-body prob- 
lems. Experimental setups where quantum states can be effi- 
ciently prepared and measured can be used to explore strong 
correlation under controlled conditions, something that typ- 
ically is not possible in the solid-state. This idea motivates 
the experimental paradigm of analogical Quantum Simula- 
tion (QS), which in recent years has been successfully im- 
plemented with ultracold atoms in optical lattices [1, 2] and 
trapped ions [3-6] . The latter have the advantage that quan- 
tum states can be prepared and measured at the single-particle 
level with great efficiency. 

A rich variety of phenomena in condensed matter physics 
is determined by the interaction of two-level systems with 
bosonic degrees of freedom. Jahn-Teller models describe 
molecules where electronic orbitals are coupled to vibrations 
such that energy is minimized by breaking some spatial sym- 
metry [7]. The many-particle extension leads to cooperative 
Jahn-Teller (cJT) models [8], which play a role in the descrip- 
tion of colossal magneto-resistance in manganites and high 
Tc- superconductivity [9, 10]. Collective effects induced by 
strong spin-boson couplings lead to intriguing phenomena, in- 
cluding effective interactions and structural phase transitions 
[7]. In this work we propose the trapped ion analogical QS 
of a one dimensional cJT model of two-level systems coupled 
to bosonic modes by a symmetric (E <g> e) interaction. Our 
proposed QS allows to explore quantum effects relevant to or- 
bital physics in solids and to use trapped ions as a testbed for 
theories to describe those systems [11]. Similar single parti- 
cle situations have been previously considered with quantum 
optical setups [12]. We study the rich phenomenology of the 
many -body case and predict that the system undergoes a quan- 
tum magnetic structural phase transition with the formation of 
a spin-boson superfluid. The latter is a quantum magnetic ver- 
sion of the classical zig-zag phase observed in ion chains [13]. 

The paper is structured as follows: (i) We introduce the 
cJT Hamiltonian. (ii) We describe the trapped ion QS of 
that model and the implementation of symmetric spin-boson 
couplings by magnetic field gradients, (iii) We calculate the 
mean-field cJT phase diagram. Within that approximation 



the formation of a spin-boson condensate is predicted, (iv) 
We compute the Gaussian quantum fluctuations (Bogoliubov 
modes) around the mean-field solution, and show that the con- 
densate is destroyed in the thermodynamical limit, but quasi- 
condensation is still possible within a range of mesoscopic 
sizes, (v) We present numerical calculations with the Den- 
sity Matrix Renormalization Group (DMRG) supporting the 
above result, (vi) Upon our knowledge on the cJT phases we 
discuss the required experimental parameters in trapped ion 
setups, (vii) Finally, in the outlook we present further inter- 
esting phenomena to be explored in either ion QS experiments 
or analogous cavity /circuit QED systems. 
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Figure 1 . Quantum magnetic structural phase transition with trapped 
ions, (a) Normal state (g < g c ). (b) Quantum magnetic zig-zag phase 
(g > gc)- For moderate chain lengths long-range zig-zag order ap- 
pears together with an antiferromagnetic spin phase. 

cJT model- We consider a chain of N spins with levels |0) 7 -, 
1 1) ; - at each site j, coupled to a chain with two boson species 
with operators ap j (/3 =x,y and ft = 1 from now on), 

H cJT =H s +H b +H c , (1) 

where H s (co z ) = (co z /2) £y Oj, and H b is the free boson term, 

H h ({A j },{tj i i})=Y g A j (fp .apj+ £ ^/(aj .a w + H.c), 
PJ PJ>i 

(2) 

where A ; - is the on-site boson energy and are boson 
hopping matrix elements. Alternatively the boson bath can 
be described in terms of normal modes, a^ n — Y,jb n japj, 
with b n j the normal mode wave-functions, such that = 
L^A^al n a P^ w i m A n the collective mode energies. The 
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last term in (1) is a Jahn-Teller E <g> e spin-boson coupling [7] 
H dg) = ^L{^K;+«L-) + ^K;+«i;)}, (3) 

that can be rewritten in terms of right and left chiral operators 

a {j = ( a tj + Myj)l& a {j = ("tj-Wyj)/^* as Allows 

H c (g) =gY< a ?( a *J + a h) + Hc - (4) 

j 

The Hamiltonian (1) is U(l)- symmetric under rotations in the 
xy plane, generated by C = Y,j( a l j a rj — a \j a ij + 

Physical implementation with trapped ions - Trapped ion 
analogical QS is ideal to implement Jahn-Teller interactions 
by controlling sideband couplings [14] with magnetic field 
gradients. Let us consider a chain of N ions of mass m, 
and charge e along the z-direction trapped by electromagnetic 
fields, see Fig. 1. Effective spins in the cJT are two internal 
electronic levels of the ions with internal energy (Oq described 
by H s (cOo). The boson bath is provided by phonons corre- 
sponding to vibrational modes of the ions in the radial (x, y) 
directions. Indeed, the position operator of ion j is 

fj = r°j z + Sr xJ x + Sr y j y. (5) 

/j are the equilibrium positions along the z-axis, and 8rp j are 
radial displacement operators, subjected to a vibrational po- 
tential that includes the radial trapping as well as the Coulomb 
interaction, 

V wib = ^mco?Y,($rpj) 2 - £ J ^ 3 (8rpj-5rpt) 2 , 
PJ j^\ r j r k\ 

(6) 

where a\ is the radial trapping frequency, and we have con- 
sidered the harmonic approximation for the Coulomb repul- 
sion. We second-quantize the vibrational Hamiltonian by 
writing Srp j = r(apj -\- j), with r = \/^/2ma\. Upon 
substitution of (5) in (6), fast rotating phonon couplings can 
be neglected in a rotating wave approximation (r.w.a.), such 
that phonon dynamics are described by H\>({g\j}, {^ u1 }), 
where G\j = G\ + ScQj is the local trapping frequency with 
ScQj = — L/(^j) e 2 / (ImCQtl/j — /f | 3 ), and the vibrational cou- 
plings become ^° ul = e 2 / (2mco t \r ( j — r^| 3 ). The r.w.a which 
leads to phonon conservation holds if G\ ^> ^° ul [15]. 

To induce symmetric Jahn-Teller couplings, we assume that 
the ion chain interacts with a time-varying magnetic field 
quadrupole 

S(x,y;t)=bf(t)(e x x-e y y), (7) 

that may be created by a conducting wire parallel to the 
ion chain [16, 17]. The magnetic coupling Hamiltonian is 
H\{t) —YtfijBj, with Jlj ji x Gj + HyC?. We assume the 
quantization axis in the z-direction such that jl x = jl y = jl. 
The spin-motion coupling can be rewritten as 

H^t) = -fibf(t) ZiSrxjaJ ~ Sr yJ oj). (8) 

j 



To control the spin-phonon couplings we use a periodic driv- 
ing f{t) = (cos v b t + cos v r t) with v b / r = coo - co z ± (cot - A). 
The latter introduce A (effective trapping frequency) and co z 
(effective spin frequency). We transform the trapped ion 
Hamiltonian to a slow-rotating frame by means of a time- 
dependent unitary transformation U(t) = e~ lHot , with //o = 
Zj{(<Do - co z )(o'/2) + (oh - A) a* pj apj}, 

H s (co z ) +H b ({A 7 -}, {^ u1 }) +H c (g) +#'(*), (9) 

where g = —flbr/Vl, and H'{t) describes fast-rotating terms 
that can be neglected as long as g <C (Qo, (Of We recover thus 
the cJT model, with renormalized local boson energies given 
now by A ; - = A + ScQt j. 

We note that (9) has the peculiarity to contain positive bo- 
son tunneling terms. This complicates the discussion that fol- 
lows, so that we implement a transformation to a staggered 
spin-boson basis, 

"/i../ ►( ' *'"/!../• tf? •( li'af. (10) 

The transformation corresponds to a 7T-rotation around the z 
axis of the odd sites of the chain. The Hamiltonian (9) is un- 
changed but the tunneling becomes now t s ^ g = (— \)i~hjf^. 
The spin ferromagnetic order in the new basis corresponds to 
a staggered order in the physical basis, in which ions alternate 
spin direction and position. 

Mean-field theory - The Hamiltonian (1) and its particular 
trapped ion realization (9) pose an intriguing quantum many- 
body problem that we approach first by a mean-field varia- 
tional ansatz. We write that ansatz in terms of chiral opera- 
tors, a\ n (£ = r,l from now on), as a product state of spins 
and displaced bosons in the collective mode basis, 

|¥MF>=0|fyfc)®^ (11) 

j 

where \BjAj) = cos(0 7 -/2)|O); + e'^J sin(0 7 -/2)|l) 7 - is a co- 
herent spin state of the spin j, and |0) e is the vacuum of each £ 
chiral mode. By minimizing the energy E = (^mfI^ I^mf), 
we arrive to a set of coupled equations for the variational pa- 
rameters 0j, (j>j, and a £in , 

a e ,n = —A-J^bnjsinOj, (12) 
ZA " j 

z 2 

(Q z tm0j = - JV^/sinG/, J L i = 2£9t^ft* 

l n 

with condition <j>j = 0/ = (j), and 0=0. The latter is an arbi- 
trary choice for a direction of spontaneous symmetry break- 
ing, such that spins are aligned in the xz plane. 

To estimate the mean-field phase diagram let us con- 
sider first periodic boundary conditions, such that b n j = 
exp(—i27lnj/N)/\^N, and couplings Jjj > 0. This last con- 
dition is met if vibrational energies A n > 0, and have a mini- 
mum at the center-of-mass mode, n = 0. We find the homoge- 
neous solution 0j = 0, with value cos0 = g 2 /g 2 if (g > g c ), 
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Figure 2. Mean-field results with N = 20 ions. Energy units such 
that co z = 1, A = 2.2. (a) Phonons per site in the chiral base. 
Black squares: homogeneous chain with g = 0.25, 0.3, 0.35 and 
tjj+i = 0.5. Blue circles: Coulomb chain with g = 0.4, 0.45, 0.5, 
and distances scaled such that at the center, f ^on = m 5 ) an d 
(c) we show results for a Coulomb (continuous line) and an homoge- 
neous (dashed line), (b) Vibrational energies of the radial collective 
modes, (c) Mean number of phonons as a function of g. The phase 
transition point at g ~ 0.25 and 0.3 is visible for the Coulomb and 
homogeneous chain, respectively. 
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Figure 3. Quantum Gaussian fluctuations F§ a n . Black, blue and 
red lines correspond to y = s (spin- waves), 1 (f-phonons) and r (r- 
phonons) respectively. We have considered homogeneous chains and 
units such that (O z = 1, A = 2.2, f/j+i = 0.5 in homogeneous chains, 
(a) N = 20 ions and different g values. As expected fluctuations ex- 
hibit a peak at the critical point. Note that F[ a ^ n } an d F{a ln } are en ~ 
hanced since spin and r-phonon operators are closer to resonance in 
(4). (b) g = 0.3 and different values of N. Fluctuations diverge with 
growing N in accordance with the Mermin- Wagner theorem [18]. 



and cos 6 = 1 if (g < g c ), where the critical coupling is g% = 
AqCO z /2. The critical amplitude is a £: o = — (gVN /2Aq) sin0, 
showing condensation into the n = mode for g > g c . Note 
that condensation above g > g c implies both ferromagnetic 
ordering of spins in the xz plane together with the collective 
displacement of the center-of-mass coordinate, thus implying 
a magnetic structural phase transition. 

To predict the phases of the cJT model in an ion chain 
we have to include finite-size corrections. We consider: (i) 
Coulomb ion chains (inhomogeneous). The separation be- 
tween ions in an harmonic linear trap, dj = r^ — r ( j_ 1 , increases 

from the center to the ends of the chain. As a result, both tj^ 
and 0\j depend on the position [15]. (ii) Homogeneous ion 
chains. This case corresponds to constant distance between 
ions. It may describe one dimensional arrays of microtraps 
and also approximates locally the description of an inhomo- 
geneous chain. Figs. 2 show results for cases (i) and (ii), with 
N = 20 ions. We present results in the staggered basis such 
that we are studying the trapped ion cJT model (9) with fff g . 
Spin-phonon condensation appears first in the center of the 
chain, something that can be used to experimentally detect 
condensation. That effect is stronger in the Coulomb chain 
case due to the inhomogeneity in the local trapping frequen- 
cies. Note that uniform displacement in the staggered basis 
corresponds to zig-zag displacement in the original phonon 
basis. 

Gaussian quantum fluctuations - Quantum fluctuations de- 
stroy long-range order in the thermodynamical limit of one- 
dimensional systems where a continuous symmetry is sponta- 
neously broken as a result of infrared divergences [18]. Still 
a mean-field theory remains a fair approximation in a range 
of mesoscopic sizes of the crystal length, thus in the most 



relevant parameter regime. To quantify fluctuations, we use 
a Gaussian approximation around the mean-field solution. 
Let us define first the fluctuation operators with respect to 
bosonic degrees of freedom, Sa e ^ n = a e ^ n — a e ^ n . Spin fluc- 
tuations are defined by means of a Holstein-Primakoff ap- 
proximation around the ferromagnetic order. For this, we de- 
fine operators Uej\0j : 0)j = |0,0) 7 -, that rotate the coherent 
spin state to a product state of spins pointing in the — z di- 
rection, erf = cos QjG z - + sin 0,of, erf = cos 0/df — sin 0/df , 
j j J j j j j j j j 

with 6j ,x = UojGj lX ul_. After the rotation, we can use the 
usual Holstein-Primakoff transformation [18] where the refer- 
ence state is taken in the a z basis, of w da I •, w Sa s /, 

J J ,J 

Gj = 2SaljSa s j — 1, valid in the limit (Sal jSa s j) ^ 1- ^i" 
nally, we substitute spin and boson operators in the cJT Hamil- 
tonian (1) and expand to second order in the fluctuation oper- 
ators, such that we get 

H G = Y,KSal n Sa £in + Y<®j Sa lj Sa sJ 

n,£ j 

+ !L^j cos Qj( Sa sJ + 8alj)(8a^ n + 8a\ n + Sa^ n + 8a] n ) 
+ 9 H^j( 5a sj - 8alj)(8a^ n - 8a\ n - (8a Vin - 8a] J), (13) 

with C0j = co z / cos 0j. Hq is diagonalized by means of a Bo- 
goliubov transformation to spin-phonon fluctuation operators 

Cm •> 

8a Yin = £ (^V^ + y «Vm) (7=s,l,r). (14) 

m=l,...,3iV 

The matrices U^ m , V^ r m , define a canonical transformation to 
a set of 3N bosonic operators c m , such that /fo = L m ®mC m c m . 
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To compute quantum fluctuations we define the vacuum 
\Q), by the condition c m \Q) =0 and define the variance per 
atom for a set of the original spin-phonon fluctuation modes, 

F {8a y , n } = JjY*( a \ 8a \,n 8a YA®)i ( 15 ) 

n 

and calculate F{g aTn y (r-phonon), F{§ ain } (1-phonon), and 

Fr§ asn \ (spin- wave) fluctuations. In Fig. 3 (a) we show that 

quantum fluctuations are smaller for larger g couplings, in 

agreement with the intuition that g ^> g c corresponds to the 

classical limit. For mesoscopic trapped ion sizes (N « 20) 

condition Fi a , \ <C 1 is satisfied, consistent with the Valid- 
ly t *y=s,r,l,nj 

ity of Hamiltonian (13). In Fig. 3 (b) we show the enhance- 
ment of quantum fluctuations a function of the system size N. 
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Figure 4. Comparison for O.P. defined in the text in a short-range 
homogeneous cJT chain with t = 0.2, and units such that A = 2, co z = 
l,andN = 20. 

Numerical calculations with the DMRG method - To vali- 
date the results obtained within a mean-field theory approach, 
we have performed DMRG numerical calculations [19]. We 
choose a maximum number of bosons per site = 20, which 
accounts for a local Hilbert space of 72 states. That local 
dimension makes our calculation computationally demand- 
ing even for short ion chains. The number of states kept in 
the reduced density matrix description (or bond-dimension) is 
D = 20, and we check that the solution has converged to an 
error of 10 -3 in all quantities presented here. 

To simplify the calculation we test the mean-field theory 
on a short range cJT model, with tj^ = —tSjj-i. Also, we 
define an order parameter in terms of long-range order, O.P. = 
Hj,kHe=i,\( a t j a e : k)/N 2 . That definition has the advantage to 
hold even without assuming spontaneous symmetry breaking. 
Fig. 4 shows that the mean-field prediction closely follows the 
quasi-exact DMRG result for a mesoscopic ion crystal with 
N = 20 ions. 

Trapped ion experimental parameters - The ability to tune 
the driving frequencies v r? b yields a large parameter control in 
the analogical QS of cJT models. Considering for 40 Ca ion 
chains with trap frequency G\ = 1 MHz, and initial energy 
splitting C0q = 20 MHz, driving frequencies can be chosen to 
get co z = 20 kHz, and A = 2.2co z . Ion separation of do « 16 
/im, would yield tj ^ =0.5(Q Z . This choice corresponds to 
the results presented in Figs. 2 and 3. We use Zeeman S\/2 



levels for \0) j and such that ji = jUb£l/2, with /Ib the 
Bohr magneton. Thus, to get to critical couplings g c « 0.2co z , 
would require magnetic field gradients b = 35 Tm _1 . Those 
values can be achieved with current trapped ion technology 
with conducting wires in planar traps. Similar spin-boson cou- 
plings may be obtained by optical forces [14]. 

To create the ground state one can follow an adiabatic 
method. We start with g = and prepare the ground state 
by laser cooling of the radial modes and pumping spins to 
\0)j. We slowly increase the coupling g. Adiabaticity holds if 
g/g <C AE, with AE the energy difference between the first 
excited and ground states. We estimate this gap from the 
energies co m of the Gaussian modes c m . For the range of 
parameters considered in Figs. 2 and 3 we find minimum 
values of AE = 0.9 and 0.4 kHz at the phase transition for 
N = 10, N = 20, respectively. Full adiabaticity would require 
thus experimental durations of the order of 10 ms, compara- 
ble to heating rates in typical trapped ion setups, something 
that could allow us to explore the effects of decoherence on 
quantum phase transitions. Even more interesting a fast ramp- 
ing of g would allow us to study the quantum Kibble-Zurek 
mechanism [20] and non-equilibrium effects in our analogi- 
cal QS. The detection of the quantum phases could be eas- 
ily performed by measuring either the phonon number or the 
spin at the end of the QS [14]. In particular, spin states in 
40 Ca ions can be read out particularly efficiently by using the 
scheme presented in [21]. The spin up state is transferred to 
a metastable D 5 / 2 state, but the spin down state remains in 
the Si/2 ground state. The illumination of the 40 Ca ion crystal 
with resonant light near 397 nm and 866 nm results in laser in- 
duced fluorescence, which is imaged on a CCD camera. Those 
ion sites which have been measured in spin up emit fluores- 
cence and appear bright while the spin down sites remain dark. 
The entire detection sequence will require 5 ms to prove the 
phase transition to the staggered order. 

Finally we discuss the effect of anisotropy in the cJT model. 
Radial trapping frequencies in each direction a\ x , G\ y may be 
tuned to be equal. Residual anisotropies in those frequencies 
can be compensated by choosing a set of sideband frequencies 
v r? b to tune radial modes to the same A in the rotating frame. 
Also, couplings g x , g y different for <7 X , o y would break the 
U(l) symmetry, but still yield interesting quantum phases, in 
particular one could study the transition from the symmetric 
case (g x = g y ) to the Ising-like g y = coupling which corre- 
sponds to a Jahn-Teller (E<g>f}) interaction [11]. 

Conclusions and Outlook.- We have studied the trapped 
ion QS of a cJT model that describes the coupling of two- 
level systems to a bath of bosonic degrees freedom. We pre- 
dict a quantum phase transition to a spin-boson condensate in 
the ground state. Experiments may allow us to explore non- 
equilibrium phenomena and the effect of decoherence. Our 
trapped ion QS involves a large number of spin-phonon de- 
grees of freedom and may challenge current numerical meth- 
ods for many-body problems. Our proposal is also relevant 
for other experimental setups such as cavity or circuit QED 
systems, where the coupling of emitters to arrays of cavities 
can be controlled to yield Jahn-Teller couplings [22, 23]. 

Acknowledgments.- We acknowledge EU projects (STREP 



5 



PICC and IP AQUTE), QUITEMAD S2009-ESP-1594, Y200200074, the Bulgarian NSF D002-90/08, DMU03/107, 
FIS2009- 10061, CAM-UCM/910758, RyC Contract and the German Science Foundation within the SFB-TRR49. 

D.R thanks Adolf o del Campo for discussions. 



[1] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 
(2008). 

[2] M. Lewenstein et al, Adv. Phys. 56, 243 (2007). 

[3] D. Porras and J. I. Cirac, Phys. Rev. Lett. 92, 207901 (2004). 

[4] A. Friedenauer et al, Nature Phys. 4, 757 (2008). 

[5] R. Islam et al, Nature Com. 2, 377 (2011). 

[6] C. Schneider, D. Porras, and T. Schaetz, arXiv: 1106.2597. 

[7] R. Englman, The Jahn-Teller Effect in Molecules and Crystals 

(Wiley, 1972). 
[8] A. J. Millis, Phys. Rev. B 53, 8434 (1996). 
[9] A. J. Millis, B. I. Shraiman, and R. Mueller, Phys. Rev. Lett. 

77, 175 (1996). 
[10] Y. Tokura and N. Nagaosa, Science 288, 462 (2000). 
[11] H. Kopel, D. Yarkony, and H. Barentzen, The Jahn-Teller 

effect, Springer series in chemical physics (Springer- Verlag, 

2009). 

[12] J. Larson, Phys. Rev. A 78, 033833 (2008); A. Bermudez, M. A. 
Martin-Delgado, and A. Luis, Phys. Rev. A 77, 063815 (2008); 



C. P. Meaney et al, Phys. Rev. A 81, 043805 (2010). 
[13] S. Fishman et al, Phys. Rev. B 77, 064111 (2008); A. Ret- 

zker et al, Phys. Rev. Lett. 101, 260504 (2008). 
[14] D. Leibfried et al, Rev. Mod. Phys. 75, 281 (2003). 
[15] D. Porras and J. I. Cirac, Phys. Rev. Lett. 93, 263602 (2004); 

X.-L. Deng, D. Porras, and J. I. Cirac, Phys. Rev. A 77, 033403 

(2008). 

[16] C. Ospelkaus et al, Nature 476, 181 (2011). 

[17] N. Timoney et al, Nature 476, 185 (2011). 

[18] A. Auerbach, Interacting electrons and quantum magnetism, 

Graduate texts in contemporary physics (Springer- Verlag, 

1994). 

[19] U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). 

[20] A. del Campo etal, Phys. Rev. Lett. 105, 075701 (2010). 

[21] C. Wunderlich etal, J. Mod. Opt. 54, 1541 (2007). 

[22] M. J. Hartmann, F. G. S. L. Brandao, and M. B. Plenio, Nature 

Phys. 2, 849 (2006); M. Leib and M. J. Hartmann, New J. Phys. 

12, 093031 (2010). 
[23] D. Porras and J. Jose Garcia-Ripoll, arXiv: 1107.2607. 



